Modeling quantum noise for efficient testing of fault-tolerant circuits 
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Understanding fault-tolerant properties of quantum circuits is important for the design of large- 
scale quantum information processors. In particular, simulating properties of encoded circuits is 
a crucial tool for investigating the relationships between the noise model, encoding scheme, and 
threshold value. For general circuits and noise models, these simulations quickly become intractable 
in the size of the encoded circuit. We introduce methods for approximating a noise process by one 
which allows for efficient Monte Carlo simulation of properties of encoded circuits. The approxima- 
tions are as close to the original process as possible without overestimating their ability to preserve 
quantum information, a key property for obtaining more honest estimates of threshold values. We 
numerically illustrate the method with various physically relevant noise models. 



Representing and transforming information using the 
principles of quantum mechanics implies that a quantum 
information processor can solve certain problems expo- 
nentially faster than any currently known classical algo- 
rithm [TH3]. Unfortunately, physically realizing such a 
processor is a difficult task as quantum systems are ex- 
tremely sensitive to environmental noise effects. Hence, 
one will likely have to take advantage of quantum error- 
correction techniques [3HZ] and fault-tolerant encodings 
of quantum information [&H10] to perform accurate large- 
scale quantum computation (QC). 

In order to verify fault-tolerant quantum computation 
is possible for a given noise model, one must prove that a 
"threshold theorem" is satisfied. The main idea of such 
a theorem is that if the error-rate on the physical opera- 
tions is below a threshold value r t h, and one can find an 
encoding that propagates errors in a controlled manner, 
then the error in the computation can be made arbitrarily 
small by concatenating the encoding. The 0-level corre- 
sponds to the physical unencoded operations, the 1-level 
corresponds to the first level of encoding and so on. 

The ability to simulate the behaviour of encoded cir- 
cuits by Monte Carlo methods can provide valuable in- 
formation such as the existence and numerical estima- 
tion of rth, as well as how errors propagate through en- 
coded operations fTTHlS] , In general, these simulations 
are inefficient on a classical computer for even moder- 
ately large systems, and so one typically makes assump- 
tions about either the types of encoded gates included 
in the circuit or the noise model, or both. In particular, 
one often restricts attention to encoded stabilizer (Clif- 
ford) circuits |17) . which is not a significant limitation 
due to the "magic state" model of QC [T5] where Clifford 
gates, ancilla magic states, and computational basis mea- 
surements suffice for universal QC. In certain cases, such 
as the seven-qubit Steane code [6], the encoded Clifford 
gates are comprised solely of 0-level Clifford gates and 
so encoded stabilizer circuits consist only of 0-level Clif- 



ford gates. By the Gottesman-Knill theorem [T7] when 
these circuits are augmented with computational input 
states and measurements they are efficiently simulatable 
on a classical computer if noise is not taken into account. 
Unfortunately when noise is taken into consideration, 
stabilizer circuit simulations are generally no longer ef- 
ficient. Hence, one often makes assumptions about the 
noise model at each faulty location, for instance, that 
the noise is described by a Pauli channel. More gener- 
ally, when the noise at each fault location is modelled 
by a mixed- Clifford channel, classical Monte Carlo sim- 
ulation is still possible by the Gottesman-Knill theorem. 
Pauli channels have various useful properties which in- 
clude a simple geometric interpretation, closure under 
composition, and diagonal x (process) matrices [H] when 
expressed in the Pauli basis. As well, Pauli channels rep- 
resent a wide class of physically realistic noise models 
such as dephasing and depolarizing processes. 

In reality, the noise at each location of the circuit will 
neither be a Pauli nor mixed-Clifford channel and so an 
important question for efficient simulation is how one can 
approximate the true noise at each location by one of 
these channels. One method is to diagonalize the noise 
in the Pauli basis by removing the off-diagonal elements 
of the ^-matrix |19j . This is an attractive method since 
in theory it can be performed experimentally via a proce- 
dure called Pauli twirling [20]. Unfortunately, in practice, 
the twirling elements will not be implemented perfectly 
and exact diagonalization can not be achieved. More- 
over, as we show later, twirling allows for the possibility 
that many states are better preserved under the twirled 
channel than the true channel, a scenario that is not ideal 
when attempting to find values of rth- 

The goal of this paper is to provide a method to ap- 
proximate a noise process A by a channel A_4 such that: 

1. A_4 is as "close" (faithful) to A as possible, 

2. A_4 provides an "honest" description of the reliabil- 
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ity of A for preserving quantum information, 

3. A,4 allows for efficient simulations of properties of 
quantum circuits. 

Point 3, coupled with the Gottesman-Knill theorem, 
motivates analyzing the specific cases of being either 
a Pauli or mixed-Clifford channel, which will be the fo- 
cus of this paper. Ideally, these A A will provide more 
realistic estimates of the threshold parameter. We em- 
phasize however that the theory developed below holds 
for completely general approximations A^. 

To quantify "close" and "honest" we will require rigor- 
ous methods for comparing quantum states and channels. 
Two standard methods for such comparisons are derived 
from the 1-norm on linear operators, || ||i, and the di- 
amond norm on linear superoperators, |j ||<> [21]. These 
comparative measures are ideal because they have a clear 
operational interpretation in terms of maximal distin- 
guishability of quantum states (operations) via POVM 
measurements . 

Our goal can now be phrased precisely in the following 
manner: 

Suppose A is a quantum channel whose complete de- 
scription is given. We want to find the Pauli, or more 
generally mixed-Clifford, channel A A that is the solution 
to the following constrained optimization problem: 

Minimize: ||A^4 — A||<> 

Subject to: for every quantum state p, 

||(A^ -^)(p)||x > ||(A-^)(p)||x. (f) 

When A_4 is a Pauli channel we denote it by A-p and 
when it is a mixed-Clifford channel it will be denoted by 
Ac- For any quantum channel A we call || (A — I) (p)\\ i 
the "input-output distinguishability" of A with respect 
to p. 

We first restrict attention to single-qubit, unital quan- 
tum channels and discuss generalizations to multi-qubit 
systems and non-unital channels later. Thus for now the 
Hilbert space of the quantum system is given by % = <C 2 
and A is assumed to map the maximally mixed state | 
to itself. The single-qubit case is particularly relevant 
when the noise affecting the circuit is highly local. The 
validity of this assumption depends on many parameters 
such as the geometric layout of the circuit, the ability 
to address specific qubits involved in a particular oper- 
ation [22], and the form of the encoding. For instance, 
since generating sets of the unitary and Clifford group 
contain only one and two-qubit gates, it may be that the 
true noise is described by a highly local model when the 
encoded operations are performed transversally as in the 
Steane code. 

Our first task is to obtain a general state-independent 
form for Eq. 0. Afterwards we will discuss restrictions 
to our cases of interest. To begin, we make use of the 



Bloch sphere representation [22] of single qubit states. 
For each state p we can associate p —> r where the 3- 
vector r lies in the unit sphere of R 3 . This vector is 
called the Bloch vector (representation) of p. For single 
qubit states, every point in the unit sphere is associated 
to a unique quantum state and the boundary (shell) of 
the unit sphere corresponds exactly to the set of pure 
states. 

Quantum channels also take a simple form in the Bloch 
sphere representation 2 1 2(1) . Any unital quantum chan- 
nel A can be uniquely represented by a real matrix M\ 
such that r M^r. This representation preserves many 
intuitive features of quantum operations, for instance a 
unitary operation IA is represented by an orthogonal (ro- 
tation) matrix Mu and Pauli channels are represented by 
diagonal matrices. 

We will exploit an extremely useful correspondence re- 
lating the 1-norm distance between quantum states p\ , p2 
and the standard 2-norm (Euclidean) distance between 
their Bloch sphere representations fi, r*2, 

1 1 Pi - Pa ||i = \\n - r 2 || 2 . (2) 

Eq. ([2]) allows us to obtain a state-independent version 
of Eq. Indeed, Eq ^ implies Eq. ([I]) is equivalent 
to 

\\r - M A r\\ 2 > \\f - M A r\\ 2 (3) 

holding for all unit vectors r. Using the theory of 
quadratic forms, Eq. ([3| is equivalent to A > B where 

A := (1 - M A f (1 - M A ) , (4) 
B :={t — M A ) T (1 - M A ) (5) 

and "T" denotes the transpose operation. In the case of 
a Pauli approximation A-p (which has a diagonal Bloch 
matrix representation), A — (1 — M-p) . Since a descrip- 
tion of A is assumed to be given, B can be computed in a 
straightforward manner [27] . Thus finding values for the 
elements of M A which give A — B > will ensure Eq. 
is satisfied. Minimizing \\A A — A||<> over these possible 
values gives the solution to our problem. 

We first look at Pauli channel approximations and 
then, using the intuition gained from these examples, 
discuss mixed Clifford channel approximations. Before 
analyzing the results in detail, let us set some notation. 
For a single qubit, the set {00, f 1, 02, 03} = {1, A, Y, Z} 
is the usual orthogonal Pauli basis for the set of 2 x 2 
complex matrices and is denoted by V\. The extension 
to the multi-qubit Pauli basis V„ is obtained by taking 
tensor products of elements of V\ . V n consists of trace- 
less (except for 1), unitary and Hermitian matrices and 
is a group when one includes phases. The Clifford group 
Clif n is defined to be the normalizer of V n and can be 
generated by the Hadamard (H), phase (S) and CNOT 
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gates applied on pairs of qubits. Note that V n is trivially 
contained in Clif„. A single-qubit mixed-Clifford chan- 
nel Ac has the form Ac(p) — ^ZiPiCipC] where the p t 
form a probability distribution, the {C{\ form a subset 
of Clifi, and p is an arbitrary mixed state input to the 
channel. In the specific case of a Pauli channel Ap we 
have A v (p) = Y?i=$Vi°iPVi- 

Numerical Results - We perform the approximation 
scheme on three types of unital channels, each of which 
represents physically relevant noise: 

A«0>) = (l-p)p + p(n p -a)p(n p -a), (6) 

3 

A {2 \p) = (l-3p)p+ P J2<npvi> (7) 



A (3 ' k) (p) = exp (-i^n k ■ a^j pexp ( 



The parameter values we choose are p = 0.01, 
n p = (sin(7r/8),0,cos(7r/8)), 9 = 0.02, and n k = 



(sin(0 fe ),O,cos(0fc)) where 9 k 



kn 



for k e {0,1,2,3,4}. 



A^'(p) represents dephasing noise about a non-Pauli 
axis, A( 2 \p) represents a depolarizing channel and the 
A^ 3 ' k \p) represent rotations about axes in the x — z plane 
starting from the z axis and ending at the x-axis. The de- 
polarizing channel was included in the analysis to verify 
that when A is itself a Pauli channel, the scheme returns 
A-p = A. The Pauli approximations A-p for each case are 
given in Table [T] via their ^-matrix xv ■ The diamond 
norm distance for each case is also given and was calcu- 
lated using the semidefinite program of Ref. [55] . 



given in Table [ill below. As expected from the constraint 



TABLE II. xo,o elements of A (1) , A (2) , A (3j) for j = 0, 1, 2. 



Channel: A« A< 2 ) At 3 '-?') 
Xo, : 0.9900 0.9700 0.9999 



of Eq. |l|, [xp]q o < Xo o m au cases. As well, since A^ 2 ^ 
is depolarizing, A^ = A^ which verifies our consistency 
check. For the unitary rotations A^ 3j \ the decrease in 
average fidelity is relatively large and is a maximum at 
k = 2. When k = 0, A^' ' is dephasing about the z- 
axis, and as the rotation axis angle approaches J, the 
approximation converges to a depolarizing channel. By 
symmetry, as the rotation axis angle goes to ?, Ap' 4 - 1 
becomes dephasing. 

It is straightforward to show there is a dephasing chan- 
nel which exactly reproduces the input-output distin- 
guishability of a unitary rotation about any axis. Indeed, 
for every rotation angle and state p whose Bloch vector 
r is at an angle a relative to z, 



| (A (3,0) _ X )(p)|| 



|sin(a)|||f|| 2 . (9) 



As well, for a dephasing channel A^' " 1 about z given by 
guishability has a very similar form: 



AJ£ (p) = (1 — p)p + pa z p<j Zl the input output distin- 



||(A^ 0) -Z)(p)|| l =2p|sin(a)|||r1 2 . 



(10) 



TABLE I. Pauli channel approximations A^ and diamond 
norm distance between Ap and A^ . 



Ap 7 - A< J 



Channel Approximation A v 

[xp]q,q \Xp\m lxv} 2 ,2 [XP] 3 ,3 

A^) 0.9860 0.0020 0.0040 0.0080 0.0152 

A< 2 ) 0.9700 0.0100 0.0100 0.0100 

A^°) 0.9900 0.0100 0.0281 

A^ 3 ' 1 ) 0.9860 0.0022 0.0040 0.0078 0.0359 

A< 3 ' 2 ) 0.9850 0.0050 0.0050 0.0050 0.0381 

A( 3 < 3 ) 0.9860 0.0078 0.0040 0.0022 0.0359 

A< 3 ' 4 ) 0.9900 0.0100 0.0281 

We give xv rather than Mp because each can be ob- 
tained from the other in a straightforward manner [27], 
the x-matrix is a standard tool in process tomogra- 
phy [T5], and any entry of the x~ matrix can be directly 
estimated via experiments [29"f[3"5] . In particular, for any 
quantum channel A, xo,o is directly related to the average 
fidelity of A, T^~i = Jtr(A(\ip)(ip\)\ip)(ip\)dip, which is 
a standard experimental figure of merit for how close an 
intended unitary operation is to the implemented opera- 
tion. The xo.o elements of A' 1 ), A^ 2 ' and the A( 3j ) are 



Hence, setting p — |sin(0/2)| implies one can exactly 
match the input-output distinguishability of A^ 3,0 ) by a 
dephasing channel. Interestingly, this dephasing channel 
is the optimal channel found by our algorithm. More gen- 
erally, a rotation about any axis n and a dephasing chan- 
nel about n have the same input-output distinguishabil- 
ity if p — | sin(0/2)| and, if the rotation is about a Pauli 
axis, this dephasing channel is a Pauli channel. 

For each j, we can use the Bloch sphere to visualize 
the difference between Ap ^ and the dephasing channel 
described above which exactly matches the input-output 
distinguishability for each p. We denote these dephasing 
channels by Ap . Fig. Ill contains plots of the defor- 
mation of the Bloch sphere in the x-z plane by Ap ^ 
(blue) and A^'^ (red) for j = 0,1,2. To make the vi- 
sualization more apparent we rotate by a larger angle, 

iverges from Ap'"^ as 9j goes to j. 
Hence, the best Pauli approximation becomes signifi- 
cantly worse than the best dephasing approximation. 
This indicates it can be useful to search over more general 
classes of channels than Pauli channels. For instance if 
we augment V\ with the Hadamard gate H then at j = 2 



= 2 sin 
Clearly A^ j) 
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FIG. 1. A^' 3) (blue), Ag j) (red), h.f' 3 > (gold) for j = 0, 1, 2. 
(Left Column): Action on the x-z plane of the Bloch sphere. 
(Right Column): Input-Output distinguishability for states p 
at angle a relative to z-axis. 

our approximation A_4 is a mixed-Clifford channel and is 
exactly Ap' 2 -*. This is because A^' 2 -* has the form 

A^ 2) (p) = (l-p)p + pHpH (11) 

where p = | sin(6*/2) | = \/0T. Thus knowing the form of 
the noise A can provide intuition as to what gates can be 
augmented to V\ to obtain a better approximation. 

Another example helps to illustrate this point. Let 
us look again at approximating A( 3,0 \ however this time 
we augment the Pauli group with the Clifford operator 
■^7r/2 = exp (— ij&z) ■ Denoting the approximation by 
Az^/2 we obtain the results in Tablelllllbelow (the results 
for using only Pauli operators are in Table]!]). The average 

TABLE III. Parameters for Approximation A^ , 2 . 

Channel Approximation Az |j Az x/2 — A' 3 ' ) || 

Xoo Xn X22 X33 

A^ ' 0.9929 0.0071 0.0151 



fidelity increased and ||A^ 7r/2 — A( 3:0 )|| o is significantly 



TABLE IV. Comparisons of ||A W - A^H* vs ||A W - A^| 



Channel 


||A«-AW|U 


||A«-4 J |L 


A« 


0.0071 


0.0152 


A( 3 .°) 


0.0020 


0.0281 


A (3-i) 


0.0020 


0.0359 


A( 3 > 2 ) 


0.0020 


0.0381 



smaller than ||A.p — A^ 3 ' ^|| o . Both of these results are 
expected since Zi provides better information for the 
approximation than Z n . 

Let us now compare some of the Pauli approxima- 
tions of Table U with those that would be obtained from 
twirling A over the Pauli group. For an arbitrary chan- 
nel A, we denote the Pauli twirl of A by A t . We an- 
alyze how A^ and A$p compare in terms of the two 
conditions we require for an optimal Pauli approxima- 
tion: first, how || — I|j compares to ||A( — for 
i = 1, (3, 0), (3, 1), (3, 2) (see Table [lV| , and second, how 
|| (A^ 3j) -l) (p)||i and || (A<?' j) -T\ (p)||i compare for 

j = 0, 1, 2 and states p at angle a relative to the z-axis 
(see right hand side of Fig. [I]) . 

We see that for all channels examined \\A^ — A^ ||<> < 
||AW — Ap || o and so the Pauli twirl is a better model of 
the original channel. However from the right-hand col- 
umn of Fig. [T] we see that, except for a small set of states, 
\\(A$ -l)(p)h > \\(A^-l)(p)\U > ||(A«-I)(p)||i 
which demonstrates that the truncation underestimates 
the error and our channel over-estimates as designed. So, 
for the price of a slightly larger diamond norm distance 
to the true error channel A, we have ensured that we do 
not underestimate the error. This can be crucial when 
considering estimates of the threshold. 

We now discuss generalizations of these results to both 
the non-unital and multi-qubit cases. The Bloch sphere 
representation of a non-unital single-qubit channel A is 
completely specified by a matrix Ma and a vector t which 
represents the non-unitality of the map, 

r^M A f+t. (12) 

One can show using a similar argument with quadratic 
forms that if A > B where 

A := (1 - M A f (1 - M A ) , (13) 
B := (1 - M A f (1 - M A ) + (||t1|| + 2||tf|| 2 ) 1 (14) 
then for every r, 

||r - M A r\\ 2 > (M A r + t) || 2 . (15) 

Hence for every state p, Eq. |l]) is satisfied where the 
vector v above is given by v = (1 — M A ) T t. 

For the multi-qubit case, it is not true in general that 
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TABLE V. Pauli channel approximation details for A^ 2q K 



Channel 


Approximation A^ 9 " 1 


HA^-A^IU 




X00 Xxx 




A( 2 <3) 


0.9900 0.0100 


0.0281 



for states p\ and p2, \\pi — P2 1| l = 11^1 — r ~2|| 2, however it is 
still likely the case that if A > B then Eq. is satisfied. 
The multi-qubit case is of significance when considering 
correlated noise models in encoded circuits. For instance, 
while many treatments of noise models in fault-tolerant 
circuits assume local, stochastic noise models, it is en- 
tirely possible that errors at certain locations can imply 
errors occur at other specific locations. It can also be 
the case that two locations always feel the same environ- 
mental influence and thus will undergo collective noise. 
Here we numerically analyze the collective unitary noise 
model A^ 2 "?) with Kraus operator exp (— iO.Olax <g> ax) 
(a two-qubit rotation about the x-axis by 6 — 0.02). The 
%-matrix of this channel is 16 x 16 but, given that it only 
contains 1 ® 1 and a x <g> a x terms, it is sparse and can be 
represented by a 2 x 2 matrix. 



0.9999 O.OlOOi 
-O.OlOOi 0.0001 



(16) 



,(2«) 



The Pauli channel approximation for this channel, A v , 

is displayed in Table [vj As expected, A^' 4 ' is symmetric 
across both qubits and interestingly gives the exact same 
results as for the single-qubit rotation about ax, Ap' 4 \ 
given in Table [TJ 

To conclude, we have provided a method for approx- 
imating quantum operations such that the approxima- 
tions are as close to the true operation as possible without 
overestimating its ability to preserve quantum informa- 
tion. We have explicitly analyzed single-qubit Pauli ap- 
proximations, discussed how one can extend the analysis 
to mixed-Clifford approximations, and have shown that 
alternative methods such as twirling the original chan- 
nel leads to highly dishonest approximations. We have 
also presented a rigorous extension to the non-unital case 
and have provided numerical evidence that our method 
likely holds in the multi-qubit case as well. These results 
are essential for simulating quantum circuits as they al- 
low for both more honest estimates of threshold values 
as well as scalable simulation of circuit properties such 
as error propagation. 
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